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Abstract 


An new eigenvalue M-linear problem arisen in the theory of metamaterials is stated 
and constructively investigated for circular non-overlapping inclusions. An asymptotic 
formula for eigenvalues is deduced when the radii of inclusions tend to zero. The nodal 
domains conjecture related to univalent eigenfunctions is posed. Demonstration of the 
conjecture allows to justify that a set of inclusions can be made neutral by surrounding 
it with an appropriate coating. 
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1 Introduction 

Local fields in fibrous composites are described by solutions of the Riemann-Hilbert and 
the M-linear problems for multiply connected domains D3 EEJ122 EH EH EH]- The physical 
properties of the components of traditional composites are expressed in terms of the positive 
constants, c.f., conductivity, permeability, permittivity etc. 

Recently, materials having negative physical constants were discovered. It concerns 
dielectric-magnetic materials displaying a negative index of refraction [UI21I91 |23l El]. Math¬ 
ematical modelling of metamaterials and neutral (invisible) inclusions were discussed in 
[121113 E] and works cited therein. In particular, the paper [12] contains a general observa¬ 
tion that any shaped inclusion with a smooth boundary can be made neutral by surrounding 
it with an appropriate coating. This result is based on the study of the eigenvalues of the 
M-linear problem for a doubly connected domain D when the spectral parameter is assigned 
only to one component of dD. Such a problem can be considered as a modification of the re¬ 
sult [29, 19] devoted to eigenvalues of the M-linear problem with the same spectral parameter 
in each component of the boundary. 

The discussed eigenvalue problem differs from the classic problem when the spectral 
parameter A enters into equation, for instance, A u + Am = 0 [fij. Our eigenvalue problem 
is similar to the Steklov problem [16] when A u = 0 in D and u = A|^ on the boundary. 
Similar mixed boundary-spectral M-linear problems were studied in ffl. 0 by reduction to 
integral equations and in [3] by variational methods. 

In the present paper, we state the general eigenvalue M-linear problem arisen in the theory 
of metamaterials and investigate it for circular non-overlapping inclusions. An asymptotic 
formula for eigenvalues is deduced when the radii of inclusions tend to zero. The nodal 
domains conjecture related to univalent eigenfunctions is posed. 

2 Statement of eigenvalues M-linear problem 

Let C = C U {oo} denote the extended complex plane. Consider n non-overlapping simply 
connected domains Dk (k — 1, 2, • • • , n) lying in the unit disk U and the multiply connected 
domain D = U\ U£ =1 (Dk U dD *,) (see Fig]l]). Let D 0 denote the exterior of the closed unit 
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disk to the extended complex plane. Let the boundary of each D k (k — 0,1, • • ■ , n) be a 
positively oriented smooth simple curve T^ including the unit circle Tq. 



Figure 1: Multiply connected domain D with circular boundaries. 

Given Holder continuous functions a k (t), b k (t ) on T fc satisfying the inequality \a k (t)\ > 
|6fc(t)| (k = 0,1,2,-•• , n). It is assumed that the winding number (index) of each a k (t) 
vanishes m- To find functions <p k {z ) analytic in D k , respectively, continuous in the closures 
of the considered domains and to find a complex constant A ^ 0 such that the following 
IR-linear conditions are fulfilled 


cp(t) = a k (t)<p k (t) + b k (t)<pk(t), teT k , k = 1,2, • • • ,n, (2.1) 

ip(t) — A oo(t)(p 0 (t) + b 0 (t)ifo(t), |f| = 1. (2.2) 

It is assumed that the unknown function tpo{z) is analytic in \z\ > 1 continuous in \z\ > 1 
and vanishes at infinity: 

^o(oo) = 0. (2.3) 

A non-zero function satisfying (12.1[) - H2.3[) is called the eigenfunction and the cor¬ 

responding constant A the eigenvalue of the problem. The function <po(z) is distinguished 
from others, since the function co(z) = (po (=), \z\ > 1, plays the key role in the theory of 
metamaterials. The univalent function u(z) determines the shapes of the inclusions (n(Tfc) 
(k — 1, 2,..., n) and of the corresponding neutral coating a;(T 0 ). 

It follows from Bojarski’s theorem [6, EJ that the eigenvalues of the problem (12.ID - (12.HD 
satisfy the inequality |A| < 1. 
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3 Functional equations 

We consider the problem fj2.ip - f)2.3p with the constant coefficients a k (t) = 1, b k (t) = —p k , 
where \p k \ < 1 (k — 1, 2,..., n) and a 0 (t) = 1, b 0 (t ) = —1. It is also assumed that T k are 
circles \t — a k \ = r k . Then, (I2.1IMI2.3D become 


( p(t) = cp k (t) - p k <Pk(t), teT k , k = 1,2, 


= *<Po(t) ~ <Po(t), |t| = 1, 

<A)(oo) = 0. 


(3.1) 

(3.2) 

(3.3) 


The problem (I3.ip - fl3.3lf can be stated in terms of harmonic functions [201 [25]. For instance, 
the condition (13. 2 p for real A np to an additive constant can be written in the form 


,, v du ., _. du 0 . . 

" = (A ~ 1) “”’ fc = (A + 1) a^’ |t| = 1 ’ 


(3.4) 


where denotes the outward normal derivative to the unit circle, u =Re</? and Uo =Rey?o- 
The R-linear problem fj3.ip ~ fl3.3p describes neutral inclusions with the contrast parameters 
p k (k = 1,2,... ,n) |ig. 

Following [5DJ El] we reduce the problem (I3.ip - fj3.3p to a system of functional equations. 
Let 


Rm) 


+ 0, r 


denote the inversion with respect to the circle \z — a k \ — r k . Introduce the function 


M*) + (z( m) ) + <P0 (=), 

jz - a k \ < r k , k — l,2,...,n, 


®( z ) '■= * <p(z) + EL 1 Pm<Pm (-* m )) + <A) ( = ) , * e D, 


k ML) + ELl PmVm [z* { 


(m) 


\z\ > 1. 


analytic in D k (k — 0,1,..., n) and D. Calculate the jump of $( 2 ) across the circle \t — a k \ = 

r k 

A k \= <F + (t) - |t - a k | = r k , 

where <F + (t) := lim z _> t zeD $ (z) , (t) := lim z _^ tzeDfe $ (z). Application of (EUD gives 

A k = 0. Similar arguments for the jump Ao of $(z) across the unit circle yield Ao = 0. 
It follows from the principle of analytic continuation that $(-) is analytic in the extended 
complex plane. Then, Liouvillc’s theorem implies that $(z) is a constant. Calculation of 
this constant as $( 00 ) and using of (13.3p yields 


$(L> = pm(p 


m V m ) 


(3.5) 


m— 1 
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The definition of $( 2 ) in \z — cik\ < r k and \z\ > 1 leads to the following system of functional 
equations 


Pk( z ) — Y2m^k P 


'Pm Pm (®m) 


+ PkPk(dk) — Po (I) , 


12 - CL k | < r k , k = 1,2,... ,n, 


\ip 0 (z) = - Em =1 P* 


Pm (^ Z ( m j ^ 'pm (&m) 


, \z\ > 1. 


Exclusion of p>o(z) from (13.61) yields the system 


Pk{z) — 'Yh m ^ k Pm 

1 s~^n — 

A 2 _.m =1 P 


Pm P m (® m ) 


+ PkPk(dk)p- 


Pm 1 + 1 Js~z ) Pm (<hn) 

- flfcl < r k , k — 1,2,... ,n. 


(3.6) 


(3.7) 


We will assume that v 7 fc(-) are analytic in |z — a^l < and continuously differentiable 
in \z — a k \ < r k due to the physical treatment of < p k (z) as complex potentials. Introduce 
the space of functions C 1 (U™ =1 .D m ) analytic in the non-connected domain U[[ )=1 D m and 
continuously differentiable in its closure with the norm 


M* = max max \ip m (z)\ + max max \p' m (z)\, 

771 1,2,... ,71 | Z CLrn\ —7*m 771. 1,2,.. .,72 | Z —7*m 

where (/>(z) = (p m (z) in | z — a m \ < r m . One can write the system (13.71) as an equation in the 
Banach space C (U^ n=i D m ) 

4> = A(j) + (3.8) 

A 

where the operators A and B are introduced in accordance with (13.71) for shortness. Equation 
(13.8p can be considered in the Hilbert space 7f 2 (U™ =1 .D m ) of functions (p(z) = <p k (z) which 
belong to the Hardy space in the disks \z — a k \ < r k with the norm [10] 

t n 1 r 2n 

V' sup — / \p m (a m + Re ie )\ 2 dd 
R < r ™ 27T 

It follows from [2D] that the operators A and B are compact in the considered spaces and 
the operator I — A is invertible where I denotes the identity operator. Then, equation (13.81) 
is equivalent to the eigenvalue problem 

A0 = (/ - A)~ 1 Bcf), (3.9) 

where the operator (/ — A)~ l B is compact in the space "H 2 (U^ =1 Z7 m ). Therefore, the eigen¬ 
value problem (13.71) can be written in the form of the eigenvalue problem (13.91) for a compact 
operator in the Hilbert space. Let 0 G "H 2 (U™ =1 .D m ) be its eigenfunction. Then, Pumping 
principle [2D] implies that 4> actually belongs to C 1 (U][ l=1 Z7 m ) (even to C°°). It is based on 
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the following arguments. For a fixed k, every function ip m (rn ^ k) is analytic in 

|z — a m | > r m and ip o (I) in \z\ < 1. The union of these domains contains the closed disk 
| z — a k \ < r k ■ ffence, the right part of (13.711 is analytic in \z — a k \ < r k . Therefore, the left 
part containing the function <Pk(z), is also analytic in \z — a k \ < r k . 

Instead of the functional equations (13.7[i we consider equations in the space C(U^ =1 .D m ) 
associated with continuous functions obtained by differentiation of (13.71) 

■0fc(~) = Yhmpk P>n ( z -a m ) 2 ^ m ( Z (m)) + 

2 , 2 \ (3.10) 

\z - a k | < r k , k = 1,2,.. . ,n, 


where ^ k {z) = p' k {z). One can see from the second equation (13.6|) that po(z) does not depend 
on <Pm(a m ). Therefore, one can first solve the system (13.101) and determine 


<Po( z ) 


i 71 2 

\ ^ Pmfm 
A “ (z — a m ) 2 

m= 1 K 



zj > 1. 


The function ipo(z) is uniquely found from (13.111) by integration 


M z ) 


/ OO 

»>'o(C)<iC, \z\ > 1. 


(3.11) 


(3.12) 


Remark 1. The eigenvalues of the Laplace operator form an increasing sequence (8]. In our 
case, the eigenvalue problem (13.91) or (13.1 0D is addressed to a compact operator. Therefore, 
the absolute values of eigenvalues decrease to zero [15]. 


4 Asymptotic solution of functional equations 

In the present section, we find asymptotic solutions of the systems (13.61) and (13.7|) when 

r 2 

r = maxfc =1) 2 ,...,n r k tends to zero. The parameters u k = are considered as values for which 
0 < z/ fc < 1 including the limit case, as r —» 0. 

Lemma 2. The eigenvalues A = A(r) satisfy the asymptotic relation 

A(r) = r 2 A 0 (r), as r —)■ 0, (4.1) 

where the function A 0 (r) is bounded as r tends to zero. 

Proof. The functions f> k {z) analytic in \z—-a k \ < r k are represented by their Taylor series 

00 / _ \ 1 

^k{z) = ^2 f 1 —— J , \z~a k | < r k , k = 1, 2,..., n. (4.2) 

i=o V r / 


6 























Here, the coefficients a \ k ^ are normalized in such a way that they are bounded as r —» 0. For 
definiteness, the eigenfunctions are supposed to be normalized as 


Mlh = ££>''h (m) l 2 = i, ( 4 . 3 ) 

m= 1 1=0 


where (p(z) = in |z — a m \ < r m . 

Using 04.21) we write equation Q3.10|) up to O considering A(r) in general form since 

its asymptotic behavior has been not known yet 


AM = r 2 Y. m ^k (iSJp 


(m) . l 

CKn + r- 




r 2 PmV-n 

A(t*) • J1TL 1 (1 CL m . 


(1 -a m z) 2 


a. 


(m) v m a[ m) z 

0 ' l-a^z 


+ 


+ °lm 


\z - a k I < r k , k = 1,2,... ,n. 
Substitute z = a k into 04.41) and reduce the order of approximation to O ( 


A k ) _ r 2 V P™ v ™ \ m ) , J±_ p m u m \m) n I _r^_ \ 

0 (a fc -a m ) 2 0 A(r) 2^im=l (1 -a^a k ) 2 0 U(r )) 5 




(m) 


k — 1, 2,..., n. 


Differentiate equations 04.4p and substitute z = a k into the result multiplied by r 


a 


(fc) 


= —2r 3 Y, 


Pm^m n ,(^) | o 7*° \ '71 prrtMrnQrn i /O 

m^k ( a k —a,m ) 3 0 ' A(r) Z^m=l (l-a^a fc ) 3 0 V \(r) 


PmVmQ'm ( w ) 

A) 3 ® 0 

k = 1, 2,..., n. 


(4.4) 


(4.5) 


(4.6) 


This procedure can be continued to get the next equations for a[ k ^ (l — 3,4,...). 

2 _ 

We now prove that cannot tend to zero as r —> 0. If it is not so, then 04. 5 p implies 

that a tends to zero as r —>■ 0. Then, equation 04.6p implies that tends to zero as 
r —» 0 and so forth ap 1 —>■ 0 for all /. This contradicts to the normalization 04.3p . 

The lemma is proved. 


It follows form Lemma [2] that the maximally possible absolute value of an eigenvalue for 
sufficiently small r can be found in the form A = r 2 /i + o(r 2 ), where p, is a non zero constant. 
Take the main terms of 04.51) and write equations up to O(r) 


(k,0) 

mo = 


EPm^m 

(1 - a^a*,) 2 ^ 

?7l=l 2 


(m,0) 

0 > 


k = 1,2, . ..,71, 


(4.7) 


where a Q fc,0) = + O(r). Introduce the matrix F whose elements have the form 


Fmk — 


Pm^n 


(1 


(4.8) 
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The eigenvalues p of the linear algebraic system (14.71) are solution of the polynomial equation 


det (pi — F) = 0, 


(4.9) 


where / stands for the identity matrix. 

If p m = p € M for any m, the matrix (14. 8 p is self-adjoint. In this case, equation (j4.9jl has 
exactly n real roots counted with multiplicity. 

The eigenfunctions can be constructed up to O(r) by (14.41) . Let p be a simple eigenvalue 
and v = (ckq 1,0 ^, ckq 2 ’°\ ..., cto"’ 0 '*) be the corresponding eigenvector of the linear algebraic 
system (ED- Then, flSUD - flOZD yield 


ipo ( z ) = -r 2 ^^^ a l ) m ’ 0 \ \ Z \>1. (4.10) 

Z O'm 

m= 1 

Example 1 (|IHj). Let n = 1 and Tx = {t e C : |f| = r} with 0 < r < 1 in the problem 
(13. ID - (13.31) . All solutions of this problem have the following form 

Vi\z) = z P , To\z) = - —, ^ (p) 0) = z p - —, A p = pr 2p , p=l,2,...,p (4.11) 

where the normalization ||</?^||^2 = 1 chosen in accordance with (14.3[) . 

The case p — 1 in Example |T] corresponds to (14.101) with n — 1, aq — 0 and = 1. 

The function, important in applications to metamaterials, 0 J\(z) = <Pq (=) = z is univalent 
in the unit disk and determines a circle neutral inclusion with an annulus coating with a 
conductivity determined by Ai = pr 2 m 

Example 2. Let n = 2, a 1 = a, a 2 = —a where a be a positive number; Ti = {t G C : 
\t — a\ = r} and T 2 = {t G C : \t + a\ = r} where a + r < 1; pi — p 2 = p. In this case, the 
system m becomes 


( 1 , 0 ) - 

= P 


( 2 , 0 ) _ 

=P 


1 (h°) | 1 (2,0) 

> Ctn I /-i , 


(1—a 2 ) 2 u 0 


(1+a 2 ) 2 u 0 


(4.12) 


1— ah 1 '® + T^af 2 ' 01 


(1+a 2 ) 2 u 0 ^ (1—a 2 ) 2 u 0 

The eigenvalues and eigenvectors of (I4.12p have the form 

_2(l + a 4 ) _ 4a 2 

Pi = P~ - nyt) V 1 = C 1 ! !); P‘2 = P-, 


(1 — a 4 ) 2 


(1 — a 4 ) 2 


, v 2 = (-1,1). 


(4.13) 


The corresponding functions oo p (z) = (=) are given by the approximate analytical formu¬ 

lae up to a multiplier 


ui(z) 


2z 


-pr 


a 2 z 2 


, a ] 2 (z) 


pr 2 


2 az 2 

1 — a 2 z 2 ’ 


One can see that the function uji(z) is univalent in the unit disks. It corresponds to the 
maximal |Ai| = r 2 \pi\. 

































5 Discussion 

The above study and examples enables us to make the following 


Conjecture. Let pk be given real numbers. Then, all eigenvalues of the problem fid. 111 - 
(1373D are real. The set of eigenvalues is countable or finite. Let |Ai| > |A 2 1 > • • •• Then the 

corresponding eigenfunctions u p (z) = iplf' 1 (=) (p — 1, 2 ,...) satisfy inequality 

wind\ z \ =1 u p {z) < p, (5.1) 


where the winding number (or index m is defined as 


wind\ z \ =1 f(z) 


1 

27TZ 



dz. 


One can see in Example Q] that 

wind\ z \ =1 Lo p (z) = wind\ z \ =1 (p^ (z) = p. 

Moreover, max/; |Afc| = |Ai| and only the corresponding eigenfunction (p^\z) is conformal in 
\z\ > 1. 

Demonstration of Conjecture for p — 1 allows to justify that any shaped inclusion with a 
smooth boundary can be made neutral by surrounding it with an appropriate coating [T2] . 

Conjecture recalls Courant’s theorem [3] outlined below. Consider for definiteness the 
Dirichlet problem u — 0 on dfl for equation A u = —Xu valid in a domain D. The set 
of eigenvalues consists of a sequence 0 < Ai < A 2 < ... (see Remark on page [6]) and the 
corresponding eigenfunctions u±,U 2 ,... constitute a complete orthonormal basis of LiifLf). 
The nodal set of a fixed u p is defined as the set {z e D : u p (z) = 0}. According to 
Courant’s theorem [3] the number of nodal domains of u p is less than or equal to p, for every 

P — 1) 2,.... 

Conjecture can be stated in terms of nodal domains of the eigenfunctions Re p ( f (z) in 
|z| > 1 of the problem (13.1H - (13.311 . Instead of (15.111 one can demand that the number of nodal 
domains of Re <Pq\z) is less than or equal to 2 p, for every p = 1,2,.... 

Let 9 e [0, 27 t) denote the argument of the complex number z. It is easily seen that the 
nodal domains of the eigenfunctions Re z~ p = \z\~ p cos p9 from Example Q] are 2 p sectors 
separated by the rays arg z = ^ where m = 0,1,..., 2p — 1. 

The general problem (I2.1H - (I2.3H and its partial case (I3.1H - (I3.3H for general curves T*, have 
been not studied yet. Even in the case of n sufficiently small circular inclusions Conjecture 
has been not proven. It is reduced to the following seemingly simple question. Let points a*, 
(k — 1, 2,..., n) lie in the open unit disk and v = (ag 1,0 ^, «o 2 ’ 0 \ • • • > a o l '°^) be eigenvectors of 
the eigenvalue problem (14.711 . For which v is the function 

n 2 _ 

«,M = yj 4 m ' 0) (5.2) 

z- a m 

m= 1 

univalent in \z\ > 1 or c v(z) = <p q (|) in \z\ < 1? Does this v correspond to the maximal \p\7 
This answer is interesting even for equal p rn and r m when the number of eigenvectors holds 
n. It solves the problem of clouds of neutral inclusions. 


9 














Acknowledgement 

The author thanks Ross C. McPhedran for stimulated discussions. 


References 

[1] A. Alu, and N. Engheta, Achieving transparency with plasmonic and metamaterial 
coatings, Physical Review E, 72, 016623, 2005. 

[2] Asatryan AA, Botten LC, Fang K, Fan S, McPhedran RC, Two-dimensional Green’s 
tensor for gyrotropic clusters composed of circular cylinders, J Opt Soc Am A Opt 
Image Sci Vis. 2014 31, 2294-303. doi: 10.1364/JOSAA.31.002294. 

[3] R. Banuelos, T. Kulczycki, I. Polterovich, B. Siudeja, Eigenvalue inequalities for mixed 
Steklov problems, arXiv:0909.5473, 2009. 

[4] Poincare-Steklov integral equations and the Ricmann monodromy problem, Functional 
Analysis and Its Applications, 2000, 34, 86-97. 

[5] A. Bogatyrev, Poincare-Steklov integral Equations and Moduli of Pants, Analysis and 
Mathematical Physics Trends in Mathematics 2009, 21-48 

[6] B. Bojarski, On generalized Hilbert boundary-value problem, Soobsch. Akad. Nauk 
Gruz. SSR, 25, No. 4, 385-390 (1960). 

[7] B. Bojarski, V. Mityushev, R-linear problem for multiply connected domains and alter¬ 
nating method of Schwarz, Journal of Mathematical Sciences, 189, 68-77, 2013. 

[8] I. Chavel, Eigenvalues in Riemannian Geometry , Academic Press, London 1984. 

[9] D J Colquitt, I S Jones, N V Movchan, A B Movchan, M Brun, R C McPhedran, 
Making waves round a structured cloak: lattices, negative refraction and fringes, Proc. 
Roy. Soc. A 2013; 469 20130218. D01:10.1098/rspa.2013.0218 

[10] P. Drygas, A functional-differential equation in a class of analytic functions and its 
application, Aequationes Math. 73 (2007) 222-232 

[11] F.D. Gakhov, Boundary Problems , Nauka, Moskow 1977 (in Russian). 

[12] P. Jarczyk, V. Mityushev, Neutral coated inclusions of finite conductivity, Proc. R. Soc. 
Lond ., A (2011), 2011 Published online doi:10.1098/rspa.2011.0230. 

[13] A.L. Kalamkarov, I.V. Andrianov, V.V. Danishevskyy, Asymptotic homogenization of 
composite materials and structures, Appl. Mech. Rev ., 62 (2009), no. 3, 030802-1- 
030802-20 

[14] M. Kerker, Invisible bodies, J. Opt. Soc. Am. 65 1975, 376-379. 

[15] A.N. Kolmogorov and S.V. Fomin, Elements of the Theory of Functions and Functional 
Analysis, Dover Publications, NY, 1999 


10 


[16] V. I. Lebedev, V. I. Agoshkov, Poincare - Steklov operators and their applications in 
analysis, Akad. Nauk SSSR, Vychisl. Tsentr, Moscow 1983 (in Russian). 

[17] G. W. Milton, S. K. Serkov, Neutral coated inclusions in conductivity and anti-plane 
elasticity, Proc. R. Soc. Lond., 457 ( 2001), 1973-1997. 

[18] V. Mityushev, Boundary value problems and functional equations with shifts in do¬ 
mains, PhD Thesis, Minsk 1984 (in Russian). 

[19] V. Mityushev, Eigenvalues of the M-lincar problems, Izvestia vuzov. Math., 11 (1992), 
35-38 (in Russian). 

[20] V. Mityushev, S. Rogosin, Constructive Methods for Linear and Nonlinear Boundary 
Value Problems for Analytic Functions. Theory and Applications , Chapman & Hall / 
CRC, Boca Raton etc 2000. 

[21] V.V. Mityushev, Riemann-Hilbert problems for multiply connected domains and circu¬ 
lar slit maps, Comp. Meth. Func. Theory, 11 (2011), 557-590. 

[22] V. Mityushev, N. Rylko, Optimal distribution of the non-overlapping conducting disks. 
Multiscale Model. Simul., 10 (2012), 180-190. 

[23] J. O’Neill, 6. Selsil, R.C. McPhedran, A.B. Movchan, and N.V. Movchan, Active cloak¬ 
ing of inclusions for flexural waves in thin clastic plates Q J Mechanics Appl Math first 
published online June 25, 2015 doi:10.1093/qjmam/hbv007. 

[24] C.G. Poulton, A.B. Movchan, N.V. Movchan, R.C. McPhedran, Analytic theory of 
defects in periodically structured clastic plates, Proc. R. Soc. A471, 2012, 1196-1216. 
DOI: 10.1098/rspa. 2011.0609 

[25] N. Rylko, Transport properties of a rectangular array of highly conducting cylinders, J. 
Engineering Math. 38, (2000) 1-12. 

[26] N. Rylko, Structure of the scalar held around unidirectional circular cylinders, Proc. R. 
Soc. A464 (2008), 391407. 

[27] N. Rylko, Effective anti-plane properties of piezoelectric fibrous composites, Acta Mech 
224, 27192734 (2013) 

[28] N. Rylko, Fractal local fields in random composites, Computers and Mathematics with 
Applications 69 (2015) 247-254. 

[29] M. Schiffer, Fredholm eigen values of multiply-connected domains, Journal d Analyse 
Mathematique, 9 (1959), 211-269. 


11 



